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^ . Abstract 

We determine the asymptotic behaviour of the number of Eulerian circuits in undirected 
simple graphs with large algebraic connectivity (the second-smallest eigenvalue of the 
Laplacian matrix). We also prove some new properties of the Laplacian matrix. 
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1. Introduction 



Let G be a simple connected graph all of whose vertices have even degree. A Eulerian 
circuit in G is a closed walk (see, for example, [2]) which uses every edge of G exactly once. 
We let Eul(G) denote the number of these up to cyclic equivalence. Our purpose in this 
paper is to estimate Eul(G) for those G having large algebraic connectivity. 

Our method is to adopt the proof given in [6 J for the case G = K n . We refer to that paper 
for the interesting history of this problem, and suggest that readers who want to understand 
our proofs carefully may find it helpful to have a copy at hand. Since the publication of [6], 
the work [3] has appeared showing that counting the number of Eulerian circuits in an 
undirected graph is complete for the class #P. Thus this problem is difficult in terms of 
complexity theory. 

Here is an outline of the paper. The asymptotic formula for Eul(K n ) and our main result 
are presented and discussed in Section 2. In Section 3 we prove some basic properties of the 
Laplacian matrix, which may be of independent interest. In Section 4 we express Eul(G) 
in terms of an n-dimensional integral using Cauchy's formula. The value of the integral is 
estimated in Sections 5 and 6, using some Lemmas proved in Section 8. We prove the main 
result in Section 7. 

2. Asymptotic estimates of the number of Eulerian circuits 

In what follows we suppose that undirected graph G has no loops and multiple edges, i.e. 

G is a simple graph. (2-1) 

We also assume that 

all vertices of G have even degrees. (2-2) 
Define the n x n matrix Q by 

-1, (vj,v k ) e EG, 
Qjk = { dj, j = k, , (2.3) 

0, otherwise 
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where n = \VG\ and dj is the degree of the vertex Vj G VG. The matrix Q = Q{G) is 
called the Laplacian matrix of the graph G. The eigenvalues Ao < Ai < . . . < A n _i of the 
matrix Q are always non-negative real numbers and Ao = 0. The eigenvalue Ai is called the 
algebraic connectivity of the graph G. (For more information about the spectral properties 
of the Laplace matrix see, for example, jl] and [7].) 

According to the Kirchhoff's Matrix- Tree-Theorem, see [5], we have that 



t(G) — — A1A2 • ■ ■ A n _i, 
n 

where t(G) denotes the number of spanning trees of the graph G. 
Let p > 1 be a real number and x G M. n . We use notation 



i/p 



For p = 00 we have the maximum norm 



max \x 



The matrix norm corresponding to the p-norm for vectors is 

IIAxIL 



We denote by \\A\ 



HS 



= SUp -r — ij-^. 

x^o \m\ p 

the Hilbert-Schmidt norm of the matrix A. 



\A\ 



HS 



\1 



■jk\ 



=1 fc=i 



If / is bounded both above and below by g asymptotically, we use the notation 

f(n) = Q klM (g(n)) , 
which implies as n — > 00, eventually 

ki\g(n)\<\f(n)\<k 2 \g(n)\. 



(2.4) 



(2.5) 



(2.6) 



[2.7) 



:2.8i 



(2.9) 



(2.10) 



When functions / and g depend not only on n, but also on other parameters we use 
notation (12. 9p meaning that condition f 1 2 . 1 j) holds uniformly for all possible values of 
The main result of the present work is the following theorem. 

Theorem 2.1. Let matrix Q be the Laplacian matrix of graph G with n vertices. Let 
conditions \2. \2. ty) hold and the algebraic connectivity Ai > an for some a > 0. Then as 
n — > 00 



Bum = e klM (2^^ vm n (I - i } 



(2.11) 



where E = \EG\, dj is the degree of the vertex Vj, t(G) denotes the number of spanning trees 
of the graph G and constants k\,k2>0 depend only on a. 
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Remark 2.1. We can replace condition Ai > an for some a > in Theorem 12.11 by the 
condition that for some a > 1/2 the degree of each vertex of G at least an. 
For the complete graph K n one can show that Ai = n and t(K n ) = n n ~ 2 . 

Theorem 2.2. (variation of Theorem 4 of [6]) As n — > oo with n odd 

Eul{K n ) = 2^rr-^n^ - lj \J (l + 0{fT 1 ^)) (2.12) 

for any e > 0. 

In fact, Theorem 12.21 is stronger than Theorem 12.11 in the case of G — K n . However, the 
asymptotic estimate of Theorem 12.11 holds for considerably broader class of graphs. 



3. Some basic properties of the Laplacian matrix 

Consider the graph G such that conditions (12. ip . (12. 2p hold. The Laplacian matrix Q has the 
eigenvector [1, 1, ... , 1] T , corresponding to the eigenvalue Ao = 0. We use notation Q = Q + J, 
where J denotes the matrix with every entry 1. Note that Q and Q have the same set of 
eigenvectors and eigenvalues, except for the eigenvalue corresponding to the eigenvector 
[1, 1, . . . , 1] T , which equals for Q and n for Q. 

Since the spectral norm is bounded above by any matrix norm we get that 

n 

An-i = ||Q|| 2 < HQH2 < HQIIi = max J2\Q jk \ = n. (3.1) 

We denote by G r the graph which arises from G by removing vertices v±, t>2, . . . , v r and 
all adjacent edges. 

Lemma 3.1. Let condition \2. 1\) holds for graph G with n vertices. Then 

n 

Ai(G) < mind,-, (3.2) 

n — 1 3 

Ai(G) > 2minc/ i - n + 2, (3.3) 

3 

\ 1 {G r )>X 1 {G)-r, (3.4) 
where Xi(G) is the algebraic connectivity of G and dj is the degree of the vertex Vj G VG. 

The proof of Lemma 13.11 can be found in pi]. 

Lemma 3.2. Let condition \2. 1\) hold and the algebraic connectivity Ai > an for some a > 0. 
Then there is a constant Cqo > depending only on a such that 

HQ-'lliHIQ-'lloo^— . (3.5) 

n 

Proof of Lemma \3.2l We consider x e W 1 such that Hx]^ = 1. For simplicity, we assume 
that \x\ \ = 1. We denote by J a the set of the indices j such that \xj\ > a/8. 
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In the case of | J a \ > an /A we have that 



a 2 



W| 2 > A /— crn/4- (3.6) 



Since the algebraic connectivity Ai > an, we get that 



n||Qx||^ > \\Qx\\ 2 > Ai||x|| 2 > era||x|| 2 > cm\j — (3.7) 



a 3 n 



In the case of | J a \ < an /A we have that 

n 

||Qx||oo > {d\ + l)|si| - ^ \ X j\ - 



i=2 

> di + 1 - \ x o\ - Yl \ x *\ - 

> di + 1 - an/4 - na/8. 



(3i 



Using again Ai > an and ( 13. 2p we get that 

n — 1 

d\ > mine?,- > an > an/2. (3-9) 

o n 

Combining (13. 7p . (13. 8p and (13. 9p we obtain that 

HQxIloo > c^nllxlloo, (3.10) 
for some constant Coo > depending only on a. * 

The following lemmas will be applied to estimate the determinant of a matrix close to 
the identity matrix I. 

Lemma 3.3. Let X be an n x n matrix such that \\X\\ 2 < 1. Then for fixed m > 2 



det(/ + X) = exp ( tr{X r ) + E m {X) \ , (3.11) 

where tr is the trace function and 

II y \\ m 

Lemma 13.31 was also formulated and proved in [6J . 

Lemma 3.4. Let the assumptions of Lemma Iff. 31 hold and all eigenvalues of X are non- 
negative real numbers. Then 

det(J -X)> exp (- . (3.13) 
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Proof of Lemma \3.4\ Using Lemma 13.31 we get that 



det(J - X) = exp 1 l ] (-±) r tr(X r )j (3.14) 

Since all eigenvalues of X are non-negative real numbers 

< tr(X r ) < tr(X) WXff 1 . (3.15) 

Hence 

det(I -X)> exp (- g ^ HXlir 1 ] > exp (-^j^) • (3-16) 



Lemma 3.5. Let the assumptions of Lemma \3."A hold. Then there is a constant c\ > 
depending only on a such that 

IdetMul < Cl ^^, (3.17) 
n 

where M\\ denotes the in — 1) x (n — 1) matrix that results from deleting the first row and 
the first column of Q = Q + J. 

Proof of Lemma 13.51 Since the algebraic connectivity Ai > an, using (I3.2p . we get the 
following estimate for the degree dk of the vertex Vk € VG. 

ii — \ 

dk > raindj > an > an/2. (3.18) 

Consider the n x n matrix X such that 

x . _ / dITT' ft fa, Vj) £ EG, (vi, v k ) $ EG and j,k^l, , g ig . 

■* k \ 0, otherwise. 1 

After performing one step of the Gaussian elimination for Q + X, we obtain that 

det(Q + X) = (di + l)detMn, (3.20) 

Since the spectral norm is bounded above by any matrix norm, we get that 

n 2 

im| 2 <||X||i<— — <-. (3.21) 
d\ + 1 a 

Since Ai > an, taking into account ( 13. 1 P - we obtain that 

\\XQ-% < \\X\\ 2 \\Q-% < A < 2 (3.22) 

a\\ a A n 

Combining Lemma [3.31 with ( I3.22p . we get that as n — > oo 

det (i + XQ-^j = exp (tr (xQ- l ^j + E 2 < exp (n-^ + 0(ra _1 )^ . (3.23) 

From (I3.20p and ( I3.23P we have that as n — > oo 

(di + 1) det Mn = det (i + XQ" 1 ) det Q < det Q exp (2/ct 2 + O^ 1 )) . (3.24) 
Since Q is positive definite, using (I3.18P in f)3.24p . we obtain (13.171) . ■ 
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Lemma 3.6. Let the assumptions of Lemma \3.S\ hold. Let G r be the graph which arises from 
G by removing vertices vi,v 2 , . . . ,v r and all adjacent edges. Then there is a constant c 2 > 
depending only on a such that for any e G (0, 1) and r <n £ 

detg(G r )> det ^ . (3.25) 
(c 2 n) 

Proof of Lemma \3.b\ We give first a proof for the case of r = 1. For our purpose it is 
convenient to use notations Q = Q{G) and Q% = Q(Gi). Note that the matrix M n that 
results from deleting the first row and the first column of Q coincides with the matrix Qi 
with the exception of the diagonal elements. In a similar way as (I3.20p we get that 

det(Q + fi + X) = (d 1 + l)detQi, (3.26) 

where X is such that 

X k ={ *TT' if ^ EG > ( Ul ' Ufc ) ^ EG and ^ 1; (3 27) 

J [ 0, otherwise. 

and Q is such diagonal matrix that 

%=ln' '*Jp> V A eEG > (3.28) 
JJ ^ (J, otherwise. 

Taking into account f 13 . 2 1 ft . we have that 

\\n + X\\ 2 <\\Q\\2 + \\X\\ 2 <l + -< -. (3.29) 

a a 



In a similar way as (|3.22|) we get that 



\\p + X)Q'% < \\Q + X\\ 2 \\Q-% < — < — (3.30) 



(3.31) 



Combining Lemma [3.31 with ( I3.30p . we get that as n — > oo 

det (l+p + X)Q~ l \ = exp Ur Up + X)Q-^ + E 2 Up + X)Q-^ > 

> exp { —n-^- — h 0(n~ l ) ] . 
\ a z n J 

From (I3.26P and ( 13 . 3 1 ft we have that as n — > oo 

(d l + 1) det Qi = det (i + p + X)Q-^ det Q > det Q exp (-3/ct 2 + O^ 1 )) . (3.32) 

Since d\ + 1 < n we get (I3.25P for the case of r = 1. 

Taking into account (13. 4 p and using r times (I3.32p we get (I3.25P for the general case. ■ 

According to (12. 4p . we have that 

t(G) = -A 1 A 2 ---A n _ 1 = (3.33) 

n n z 

where t(G) denotes the number of spanning trees of the graph G. 
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Lemma 3.7. Let the assumptions of Lemma \3.2\ hold. Then for some C3 > depending only 
on a the number of spanning trees of the graph G with maximum degree greater than d is 
less then C3 det Q/d\ for all d > 0. 

Proof of Lemma \3. r A According to Lemma 5 of [6j the number of labelled trees on n vertices 
with first vertex having degree greater than d is less than 2n n ~ 2 /d\ for all d > 0. We have 
that 

detQ > A™ > (an) n . (3.34) 

To complete proof it remains to note that the number of spanning trees with maximum degree 
greater than d in G does not exceed the number of such spanning trees in the complete graph 
with n vertices. ■ 



Consider a spanning tree T of the graph G. We denote by Gt the graph which arises 
from G by removing all edges of the tree T. 

Lemma 3.8. Let the assumptions of Lemma \3.2\ hold. LetT be the spanning tree of G with 
the maximum degree at most an/ '4. Then the algebraic connectivity Xi(Gt) > an/2 and 

det Q(G T ) > c 4 det Q(G) (3.35) 

for some C4 > depending only on a. 

Proof of Lemma \3.8l Note that 

Q{G T ) = Q(G)-Q(T). (3.36) 

Since the spectral norm is bounded above by any matrix norm and the maximum degree of 
vertex of T at most an/ '4 we get that 

liQ(T)|| 2 <||Q(T)||i< an/2. (3.37) 

Therefore, since the algebraic connectivity Xi(G) > an, we have that 

Ai(G r ) > \i(G) - \\Q(T)\\ 2 > an/2 (3.38) 

and 

det Q{Gt) = det Q(G) det(I — X), (3.39) 

where X = Q{T)Q{G)^ 1 . Note that Q(G) is the matrix of positive definite quadratic form 
and Q(T) is the matrix of quadratic form with non- negative eigenvalues. Considering the 
basis in which both matrices are diagonal, we have that 

tr{Q{T)Q{G)- x ) < tr(Q(T))\\Q(G)-% (3.40) 

and 

all eigenvalues of Q(T)Q(G') _1 are non-negative. (3-41) 
Using again the fact that the algebraic connectivity X\(G) > an and ( I3.37P we get that 

\\XW2 < \\Q(T)\\ 2 \\Q(G)-% < ^— = \. (3.42) 

2 an 2 
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Since T is the spanning tree tr(Q(T)) = 2(n — 1). Using (I3.40p . we get that 

tr(X) < tr(Q(T))\\Q(G)-% < In— = -. (3.43) 

an a 

To complete the proof it remains to combine Lemma [3.41 with (13.39}) . (I3.42[) and (I3.43H . ■ 



Lemma 3.9. Let a > and the assumptions of Lemma \3.S\ hold. Then for any set A C VG 
such, that \A\ > an, there is a function h : VG — > No, having following properties: 

h(v) = 0, ifveA, h(v) < H, for any v G VG, (3.44) 

{w G VG | (w,v) G EG and h(w) < h(v)} > an, if v £ A, (3.45) 
where constants H, a > depend only on a and a. 

Proof of Lemma \3.iA At first, we construct the set A% = {uG VG \ h(v) = 1}, having 
property (13.451) . 

If \A\ > n — an/ 4, then let A\ — {v G VG \ v ^ A} . Taking into account (13. 18 j) . we get 
that property (13.451) hold for a = a/A. In this case H = 1. 
For \A\ < n — an/A define x G M. n such that 

x = { l-\A\/n, Vj G A, 

Xl \ -\A\/n, Vj <£A. [6M) 

Since %\ + x<i + . . . + x n = 

x T gx = x T gx> Ai||x||l > Ai|A| (— -i^lV >anan(a/A) 2 =°^-. (3.47) 



n J 16 
On the other hand, 

x T Qx= J2 ( x J- x k) 2 , (3-48) 

(vj,v h )eEG 

which is equal to the number of edges (v , w) G EG, where v G A, w ^ A. We denote v4x the 
set of vertices w ^ A, having at least an adjacent vertices in A, where a = -^aa 3 . 

x T Qx < n\A x \ +an\VG\. (3.49) 

Combining ( 13.47P and (I3.49p . we get that \Ai\ > an. 

We make further construction of the function h inductively, using for the k-th step the 
set A^ = A U A\ U . . . U A\~. The number of steps does not exceed 1/a as \A^\ > an for each 
step, perhaps with the exception of the last one. ■ 
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4. The result expressed as an integral 

The reasoning below is similar to the arguments of Section 2 of [6J 

An Eulerian orientation of G is an orientation of its edges with the property that for 
every vertex both the in-degree and the out-degree are equal. Any Eulerian circuit induces 
an Eulerian orientation by orienting each edge in accordance with its direction of traversal. 

A directed tree with root v is a connected directed graph T such that v G VT has out- 
degree zero, and each other vertex has out-degree one. Thus, T is a tree which has each edge 
oriented towards v. 

Let D be a directed graph with n vertices, and let v G VD. A directed spanning tree of 
D with root v is a spanning subgraph of D which is a directed tree with root v. 

The following famous theorem, sometimes called the BEST Theorem, is due to de Bruijn, 
van Aardenne-Ehrenfest, Smith, and Tutte [HE]. 

Theorem 4.1. Let D be a directed graph with vertices v i, V2, ■ ■ ■ , v n . Suppose that there are 
numbers d±, o^, • • • , d n such that, for every vertex v r , both the in-degree and the out-degree of 
v r are equal to d r . Let t r = t r (D) be the number of directed spanning trees of D rooted at v r . 
Then t r is independent of r, and 

n 

Eul(D) = t r Y\_(dj — 1)!. (4.1) 

3=1 

Consider the undirected graph G with n vertices such that conditions (12. ip . (12.21) hold. 
Note that for every spanning tree T of the graph G and any vertex v r G VG there is only 
one orientation of the edges of T such that we obtain a directed tree with root v r . We denote 

n 

by % the set of directed trees with root v r obtained in such a way. For T G |J % denote by 

r=l 

EO(T) the number of Eulerian orientations of G that the corresponding graphs contain T. 

From Theorem 14. II in the case of a graph D corresponding to Eulerian orientation of the 
graph G we find that 

Eul(D)=t r (D)H M-l !, (4.2) 

3=1 ^ ' 

where dj is the degree of the vertex Vj G VG. Let denote by SO the set of all graphs 
corresponding to Eulerian orientations of the graph G. Grouping Eulerian circuits according 
to the induced orientations, we obtain that 

Eul{G) = Eul(D) = n(f-l)<£ tr( D ) ( 4 -3) 

DeEO j=l ^ ' Deeo 

for any fixed natural number r < n. 

Regrouping the terms of the final summation according to the directed subtrees rooted 
at v r , we find that 



Eul(G) = f[(^-l)\J2EO(n 

o—1 \ / TcT. 



(4.4) 

3=1 v ~ 7 Te% 

For n > 1 and R > we use notation U n (R) = {(xi, X2, ■ ■ ■ , x n ) \ \xi\ < R for all i}. The 
value of EO(T) is the constant term in 

n (xr^k+xk-^j) n i^;';i„ r ( 4 - 5 ) 

{v h v k )&EG ( Vj ,v k )eET V 3 * K JJ 
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which we can extract via Cauchy's Theorem using the unit circle as a contour for each 
variable. Making the substitution Xj = e ldj for each j, we find that 

n j 

Eul{G) = J] (- " 1)! 2 |£Ghn+1 7T-"5, (4.6) 

where 

S= JJ cosA jfc ^ Yl (l + *tanA#) d0, (4.7) 

uJ-n/2) ivj,v k )&EG TeTr(vj,v k )€ET 

having put Ajk = 9j — 6k and using the fact that the integrand is unchanged by the 
substitutions 9j — >■ 9j + it if condition ( 12. 2 p holds. 

We approach the integral by first estimating it in the region that would turn out to be 
the asymptotically significant one. Then we bound the integral over the remaining regions 
and show that it is vanishingly small in comparison with the significant part. 



5. The dominant part of the integral 

In what follows, we fix some small constant e > 0. Define 

Vq = {6 E U n (n/2) : \9j - 6\ (modn) < n~ 1/2+£ , where 9 = °_l±i__±^^ ( 51 ) 

and let So denote the contribution to S of 6 G Vq. Since the integrand is invariant under 
uniform translation of all the 9/s mod tt, we can fix 9 = and multiply it by the ratio of its 
range n to the length n~ 1 / 2 of the vector ^[1, 1, ... , 1] T . Thus we have that 



S = -nn 1 ' 2 



I j ] cosA jfc ^ II (1 +«tanA ifc ) rfL, (5.2) 



LnV (vj,v k )eEG T&Tr {vj,v k )^ET 



where L denotes the orthogonal complement to the vector [1, 1, ... , 1] T . 

The sum over % in the integrand of (15. 2p can be expressed as a determinant, according 
to the following theorem of [9J. 

Theorem 5.1. Let Wjk (1 < j, k < n, j ^ k) be arbitrary. Define the n x n matrix A by 

the sum being over 1 < r < n with r / j. For any r with 1 < r < n, let M r denote the 
principal minor of A formed by removing row r and column r. Then 

det M, = ] [ w jk , (5.4) 

T ( Vj ,v k )eET 

where the sum is over all directed trees T with VT = {vi, t>2, . . . , v n } and root v r . 
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Lemma 5.1. Let the assumptions of Theorem \2.1\ hold. Let Q = Q + J , where J denotes 
the matrix with every entry 1. Then for 6 G Vq as n — )■ oo 

EE II (l + itauA ifc ) = detQ^l + O^- 1 / 2 ^)], (5.5) 

where at denotes the vector composed of the diagonal elements of the matrix Q^ 1 , A denotes 
the diagonal matrix whose diagonal elements are equal to the components of the vector Q6. 

Proof of Lemma \5.1\ Define the n x n matrix B by 

-tanAjfc, for (Vj,Vk) G EG, 

b ' L \ l:(vj,Vi)eEG W 

otherwise . 



Using Theorem 15. II with the matrix A = Q + iB, we get that 

n n 

EE II (l + *tanA ifc ) = X>r> ( 5 - 7 ) 

r=l TeTr (vj,v k )EET r=l 

where M r denotes the principal minor of A formed by removing row r and column r. Since 
the vector [1, 1, ... , 1] T is the common eigenvector of the matrices Q and B, corresponding 
to the eigenvalue 0, we find that 

f. det(Q + «) 

/ — ' n 

r=l 

Note that for G V 

A jk = (6 3 -6)- (6 k -&)< 2n- 1 ' 2+£ . (5.9) 
Since the spectral norm is bounded above by any matrix norm we get that 

n 

\\B\\ 2 < ||£||i = m&x^2\B jk \ = 0(n 1/2+£ ). (5.10) 
3 k=i 

Let $ = BQ^ 1 . Since the algebraic connectivity Ai > an, we get that 

11$ || 2 < 1 1^1 12I |<0 1 1 Is < ^-\\B\\ 2 = 0(n- x / 2+£ ). (5.11) 
Using Lemma (13. 3D with the matrix we find that as n — > 00 

det(J + i$) = exp f*r(i$) + + 0(n- 1/2+3e )^) . (5.12) 

Let 

B = B skew + B diag , (5.13) 

where B skew is the skew-symmetric matrix and B^ ag is the diagonal matrix. Since Q is the 
symmetric matrix 

tr(B skew Q- 1 ) = 0. (5.14) 
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Using (15. 9p . note that as n — > oo 



\B 



diag 



A I 



0{n 



(5.15) 



where A denotes the diagonal matrix whose diagonal elements are equal to the components 
of the vector QO. Since the algebraic connectivity Ai > an, we get that as n — > oo 



-l/2+3e\ 



(5.16) 



tr [{B diag - A)Q-^ \ < n\\B diag - A\\ 2 \\Q'% = 0{i 

Using f)5.14p and (15.16)) . we obtain that as n — > oo 

tr($) = HB^Q- 1 ) = triAQ- 1 ) + C^n" 1 / 2 ^) = 8 T Qa + C^rT 1 ^), (5.17) 

where a. denotes the vector composed of the diagonal elements of the matrix Q~ x . 
Using the property of the trace function tr(XY) = tr(YX), we have that 

tr ($ 2 ) = tr {{B skew Q- 1 ) 2 ^ + tr ({B diag Q- 1 ) 2 ^ +2tr (b skew Q~ l B diag Q- 1 ^ 
Since B skew is the skew-symmetric matrix and Q~ l B diag Q~ l is the symmetric matrix 



t r {B s kewQ 1 B diag Q 1 j — 0. 



One can show that 



Therefore we get that 



|2 



tr (X 2 ) < \ \X\ 
\XY\\ HS < \ \X\\ H s\\Y\\ 2 . 



tr ( (B skew Q 



-1\2 



I \ B s kewQ | \ hS 



Since the algebraic connectivity X 1 > an, using (15. 9p . we obtain that as n — > oo 

-l/2+e\ 

11^1 | J -'s/ceiu| \ ns — \ 1 1 J -'snew\ \tia — )■ 

Using ( 15. 9ft and ( 15. 15lt . we get that 



\B s kewQ 1 \\HS < \\Q 1 \ b| \B s kew\ \hS < T - 1 1 B skew \ \ }jg — 0(n 



(5.18) 

(5.19) 

(5.20) 
(5.21) 

(5.22) 



tr [{B diag - A)Q 1 B diag Q ^ < n-^\\(B diag - A)\\ 2 \\B diag \\ 2 = 0{n 



and 



tr ({B diag - A)Q-\B dmg - A)Q-') \ < n^\\{B diag - A)\\ 2 2 = 0{n 



-2+6e\ 



Thus 



tr ( (B diag Q 



-1\2 



A? 

tr [ (AQ- 1 ) 2 



0(n 



-l+4e> 



Combining (IBTTSl) . ( Kim . 1)521) . (ET221) and (15^51) . we obtain that 

tr ($ 2 ) = tr ((AQ- 1 ) 2 ) + 0(rT 1+4e ). 



(5.23) 

(5.24) 
(5.25) 

(5.26) 
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Using (l5TTTj) and flQUD in we get that 



det(7 + z$) =exp i T <2a + 



tr (AQ 



-1\2 



+ 0(n 



(5.27) 



Combining (j577j) . fl5T8|) and fl5T27j) . we obtain fl5T5|) . 

We denote by P(#) the orthogonal projection onto the space L. 



Lemma 5.2. Let the assumptions of Theorem \2.1\ hold. Let Q = Q + J, where J denotes the 
matrix with every entry 1. For positive constants a,b,c,di,d,2 let sequence of vectors {oc n } 
and sequence of differentiable functions {R n (0)} be such that 



|«n|U < c/n, 
T Q0 



\Rn(9)\<d v 



n 



R n {0) = Rn(P(0)) 



and for 0E £4(^ 1/2+E ) 



For n > 2 define 



dR n {0) 



00 



<d,n- l ' 2+e . 



J n = j exp U0 T Qa n -a ^ A%-b ^ A% + R n (0)\dL. 

LnV \ (v h v k )eEG (v h v k )eEG J 



Lnv 

Then as n — > oo 



where constants ki, k 2 > depend only on a, b, c, d\, d 2 and a. 

Lemma [5.21 is proved in Section 8. 
Lemma 5.3. Let the assumptions of Theorem \2.1\ hold. Then as n — )■ oo 



n— 1 n+1 



S = 6fci,*a 2 a 7r » n VdetQ , 



(5.28) 

(5.29) 
(5.30) 

(5.31) 



(5.32) 



(5.33) 



(5.34) 



where constants k%, > depend only on a. 

Proof of Lemma \5. 3[ Using formula (15. 2p with r = 1, 2 . . . , n and summing, we obtain that 



n „ 

cos Aj k ]I + itanA 3 fc) <iL 

r=l ^ fe,» fe )eBG Te%{v jt v k )eET 



/n 
Yl cosA jfc ^^ J~J (1 + z tan A jk ) dL. 



(5.35) 



Lnv (fj.«fc)eEG 



r=l TS7; (vj,v k )£ET 
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By Taylor's theorem we have that for 6 G Vq 

J] cosA ife = exp(-l Yl A %'Y2 S Ajk + Oin- 1 **)]. (5.36) 

(vj,v k )eEG \ {v 3 ,v k )(LEG (vj,v k )eEG J 

Combining f!5.35j) with (I5.36P and Lemma 15. 1[ we obtain that as n — > oo 

LnVo V (vj,v k )eEG {v 3 ,v k )&EG j 

where ex. denotes the vector composed of the diagonal elements of the matrix Q" 1 and 

R{0) = tr{KQ~ 1 ) 2 , (5.38) 

where A is the diagonal matrix whose diagonal elements are equal to the components of the 
vector QO. Since vector [1, 1 ... , 1] T is the eigenvector of Q, corresponding to eigenvalue 0, 
we have that 

QO = QP(0). (5.39) 

Thus 

R(9) = R(P(0)). (5.40) 
Note that Lemma 13.51 implies as n — > oo 

||«||oo < ci/n, (5.41) 

where c\ = c\{a) > 0. Since the algebraic connectivity Ai > an, using (15.20J1 . we get that 



(5.42) 



\R(0)\<\\AQ-^ s <\\M\h\\Q-% = 

= \\Q0\\l \\Q-% < \\Q\\i \\Q-% Mil < 

n l X 2 i Af Ai a n 

Note that for 6 G U n (^n~ l l 2+e ) and for some 1 < k < n 

n 

iiaii 2 <Ei a ^i = °K /2+£ )- ( 5 - 43 ) 

For 1 < /c < n we denote by (Q~ l ) k the fc-th column of the matrix Q^ 1 . Using again Ai > an, 
we get that 

1 = WQiQ^hWz > AxIKQ- 1 )*!^ > an\\(Q-%\\ 2 . (5.44) 

Note that 

— = 2tr(—Q- 1 AQ- 1 ) = 2d k (Q~ 1 AQ- 1 ) kk + 2tr{AQ~ l AQ- 1 ), (5.45) 
ov k ou k 
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where (Q~ 1 AQ~ 1 ) kk denotes (k, k)-th element of the matrix Q~ 1 AQ~ 1 and the matrix A is 

such that for any 1 < j < n the diagonal element \Ajj\ < 1 . Since the algebraic connectivity 
Ai > cm, using (I5.43P ( I5.44p . we get that as n — > oo 

ItriAQ-'AQ- 1 )] < nWAQ^AQ- 1 ^ < n\\A\\ 2 ||A|| 2 = 0(n~ 1/2+£ ) (5.46) 

and 

2d k (Q- 1 AQ- 1 ) kk = 24A fcfe ||(<T 1 ) fc || 2 , = 0{n- l ' 2+£ ). (5.47) 
Thus for some d > 0, depending only on a 

dRJO) 



< dn~ 1/2+£ . (5.48) 

of) 

Combining fl5T37j) . (00]) . (OT]l . fl5Tl2|) . f l5^8|) and using Lemma [El we obtain fOD . ■ 
6. The insignificant parts of the integral 

In this section we prove that So contributes almost all of S, even though it involves only a 
tiny part of the region of integration, compare with Section 4 of [6]. We continue to use the 
same value of e as in the previous section. 

Let assumptions of Theorem 2.1 hold. Define E'T = {(v j, v k ), {v k , v,j) \ (vj,v k ) G ET}. 
We express the integrand of (14.71) as 

f ^= e n ^( t >0)> (e- 1 ) 

TgT r jkeEG 

where 

{cosA ifc (l + itanA jk ), (vj,v k ) G ET, 

cosA ifc (l - ^ tan A jk ), {v k ,Vj) G ET, (6.2) 

cosAjfc, otherwise. 

Note that \fj k (T, 0)\ < 1 for all values of the parameters. One can show that 

1 9 
|cos(a;)| < exp(— -x 2 ) for \x\ < j^ 71 "- (6-3) 

Divide the interval [— |7r, |7r] mod n into 32 equal intervals Hq, . . . , H^i such that Hq = 
[— ^7r, ^71"]. For each j, define the region Wj C U n (n/2) as the set of points having at least 
-^n coordinates in Hj. Clearly, the Wj's cover U n (n/2) and also each Wj can be mapped to 
Wq by a uniform translation of the 6j mod 7r. This mapping preserves the integrand of (14. 7p 
and also maps Vq to itself, so we have that 

j \F(9)\dO <32Z, (6.4) 

C7„(7r/2)-V 

where 

Z= J \F{6)\dO. (6.5) 

w -v 
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We proceed by defining integrals Si, . . . , S4 in such a way that Z is obviously bounded 
by their sum. We then show that Sj = o(Sq) for j = 1, 2, 3, 4 separately. Write 

F(0)=F a (0) + F b (0), (6.6) 

where F a (9) and F b (6) are defined by restricting the sum to trees with maximum degree 
greater than an/ '4 and no more than an/ '4, respectively. Also define regions Vi and V2 as 
follows. 

Vi = {0 G Wq I I #7 1 > — vr for fewer than n £ values of j}, 



32 

V2 — {0 G Vi I > — 7r for at least one value of j}. 



(6.7) 



16 

Then our four integrals can be defined as 



Si= j \F(0)\dd, 

W0-V1 

S 2 = I \FJ0)\d0, 



Vi 

Ss = J \F b {0)\d0, 

S 4 = J \F b (0)\d0. 

Vi-Va-vb 



(6.8) 



We begin with Si. Let h be the function from Lemma 13791 for the set A = {vj \ \9j\ < ^vr}. 
We denote l m in such natural number that inequality 



J 3 



>l w (l + //fT) (6.9) 



holds for at least n e /H indices of the set {j \ h(vj) = /}. Existence of l min follows from the 
definition of the region Vi. If 6j and 9k are such that 

\0j\ > ^-Ml + Imin/H) and \9 k \ < i-7r(l + (l mm - l)/H) (6.10) 
64 64 

or vice versa, but (vj,Vk) E'T, we have that \fj k {T,0)\ < cos(^7r/iJ). This includes at 
least {an — n 5 )^ — n edges (vj, Vk) G EG. Using (13. ip and ( I3.33p . we get that as n — > 00 



5i < t(G)n n [cos —) =0^exp(-cn 1+£ )j2— vr^-n-ydetQ (6.11) 

for some constant c > depending only on a. 

To bound S2, we first note from Lemma \3. 71 that the number of trees with maximum de- 
gree greater than an/4 is less than C3 det Q/ (an/4)\. Using (16. 3p . we see that 

\f ]k (T,0)\<exp(-l-A%) (6.12) 
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except for at most n 2e pairs (J, k) with \Ajk\ > j^tt and fewer than n pairs in E'T. In those 

.lA 2 



16 

excluded cases the value exp(— | A 2 fc ) may be high by a factor exp(|7r 2 ). Hence, we have that 



s ^W§. ew ^ n+n2e) ) j exp i-^ >^ ^ d °- (6 - i3) 

U n (n/2) 

Lemma 6.1. Let the assumptions of Lemma \3.2 hold. Then 

-1 \ o— S±l 

•- £ A% }d0<— ] ==r-. (6.14) 

(v jt v k )eEG J ydetQ 

Lemma [67TI is proved in Section 8. Combining ( 16. 13[) and (I6.14p . we obtain that as n — > oo 



\ (v jt v h )eEG ) 



exp 

U n (ir/2) 



S 2 = 0{n- cn )T^Ti^n- x \JfeiQ (6.15) 

for some constant c > depending only on a. 

We denote by Gt the graph which arises from G by removing all edges of the tree T. 
Let Gt,o be the graph resulting from from Gt by removing vertices, corresponding to those 
values of j for which \6j\ > j^tt. 

For 1 < r < n £ let S^{r) denote the contribution to S3 of those 9 G Vi such that \9j\ > j^tc 
for exactly r values of j. If \9j\ < -^tt and \9k\ > ^7r or vice versa, we have that 

|/ ifc (T,0)|<cosQ^ (6.16) 

unless (vj,Vk) G -E'T. This includes at least r[an/2 — an/4 — n £ ) pairs (j,k), because the 
degree of any vertex of the graph G is at least an/2, see (13. 2p . For pairs (j,k) such that 
\0j\, \&k\ < ^vr, but (Vj,v k ) <£ E'T, we use (gTTJD - We put 0' = (0 1} . . . , # n _ r ). Then, 

S 3 (r)<^(cos^) rW4 " n£) ^^ / exp(-i £ A| fe ] <W, (6.17) 

where the first sum is over choices of those values of j for which \9j\ > ^7r and the second 
sum is over trees with maximum degree an/ '4. Using Lemma [3.81 and then Lemma [3.11 and 
Lemma [3.61 for the graph Gt, we obtain that 

\i(G T ,e) > o-n/2-n e (6.18) 

and 

detQ(G T , )>-^, (6.19) 
(c 5 n) 

where C5 = c§(a) > and Q = Q(G). According to Lemma [6.11 we have that 

exp -I £ A] t Ly< 2 ' \ * " . (6.20) 
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Combining (I6.17P with ( I6.19p . (I6.20p and allowing n r for the choice of those values of j for 
which \9j\ > j^n, we obtain that 



S 3 (r)<2 2 7T 2 n ^i(oos-j ,- = (6.21) 

'detQ 



and, using fl3.33[) . we can calculate that 



S 3 = J2s 3 {r) = 0{c- n )2^ 1 Ti^n- 1 \JdetQ (6.22 



r=l 



for some constant c > 1 depending only on o. 

Since < |-7r for G Vi — V2 — Vq and the integrand is invariant under uniform 
translation of all the Oj's mod it, we can fix 9 = and multiply it by the ratio of its range 7r 
to the length n~ 1//2 of the vector -[1, 1, . . . , 1] T . Thus we get that 

S 4 <nn 1/2 J \F b {0)\dL, (6.23) 

LnU n (ir /8)-V 

where L denotes the orthogonal complement to the vector [1, 1, . . . , 1] T . In a similar way as 
flBTTTD we find that 

S A <nn^J2 J ex P|4 E A %) dL i ( 6 - 24 ) 

T inC/n(7r/8)-Vb V (^^)&EG T I 

where the first sum is over trees with maximum degree an/ 4. 
Lemma 6.2. Let the assumptions of Lemma \3.S\ hold. Then as n — )■ 00 



/ exp ~2 E A % ) dL = O {exp(-cn^)) T =r~ (6-25) 



L-C/„(n- 1 /2+ £ ) 

/or some c > depending only on a. 

Lemma IS~2l is proved in Section 8. Using Lemma [3781 and combining (I6.25p . (I6.24p and (I3.33p . 
we obtain that asn-^oo 

~ n— 1 n + 1 

2 2 7T 2 U „ / , o_,x 1-1 n+1 



S 4 = (exp(-cn 2£ )) t(G) ==- = O (exp(-cn 2£ )) 2~ tt ^n~y det Q (6.26) 

'detQ 



for some c > depending only on o. Combining (16. lip . (16. 15ft . (16. 22ft . (I6.26P and Lemma 
15. 3[ we obtain the desired result. 

Lemma 6.3. Let the assumptions of Theorem \2.1\ hold. Then as n — > 00 

S =(1 + (exp(-cn 2e ))) S (6.27) 
for some c > depending only on a. 
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7. Proof of Lemma 13.21 

According to (l4~6|) and ( 14771) 

n , 

Eul{G) = TT ( J - l)\2^- n+1 n- n S, (7.1) 

where 

S= n cosA ifc$^ ] [ (1 + ztanA^) d0. (7.2) 

Combining Lemma 15.31 and Lemma 16.31 we get that as n — >■ oo 



n — 1 n + 1 



5 = Q klM I 2— TT— n-^detQ 1 , (7.3) 
where constants fci, &2 > depend only on a. Taking into account ( I3.33P we obtain (12. lip . ■ 



If for some a > 1/2 the degree of each vertex of the graph G at least an, we can use (13. 3 j) 
and get that 

Ai(G) >2nund i -n + 2> (2<r - l)n. (7.4) 



.7 



8. Proofs of Lemma 15.21 Lemma 16.11 and Lemma 16.21 

Let assumptions of Lemma 13.21 hold. We define 

</>=<K0) = (0l(0),...0n(0)) = Q0. (8.1) 

We continue to use notation P(6) for the orthogonal projection onto the space L, where L 
is the orthogonal complement to the vector [1,1,..., 1] T . For any a > we have that 

J e -ae T Q9 de = 7r n/2 a -n/2^ q ^ 



and 



J e- a9T ® 9 dL = J e- a0TQe dL = n^a^n 1 / 2 / yftetQ. (8.3) 

L L 

Proof of Lemma \6.1\ Note that 

&% = T Q6. (8.4) 

(vj,v k )£EG 

Since diagonal of U n (Tc/2) is equal to 7m 1 / 2 and QO = QP{0) we have that 

exp j -I Y, A % ) de < ™ 1/2 [ e-^ 0TQe dL. (8.5) 



U n (w/2) 

Using (IQjl . we obtain ( 167141) 
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Note that for some g%{9) — gx(02, ■ ■ ■ , 9 n ) 



T Q0=P°¥-+ gi (e). (8.6) 



Using (13. 9p . we get that as n — > oo 

+oo 

J e- aeT WdO = j ■■ j e- a9 ^>~^ I j e~ a ^d9 1 J d0 2 ...d6 n 

R n -oo -oo Voo / 

(l + 0(exp(-m 2£ ))) J e- aeT ® e dO 



-oo +oo / +oo 

,",2 



•7) 



l (6)\<hc£n 1 / a +' 



for some c > depending only on a and a, where is the constant of Lemma I3T21 Combining 
similar expressions for <fti, <p2, ■ ■ ■ <fin, we obtain that as n — > oo 

J e - aeT Q e dO= (l + 0(exp(-cn 2£ ))) J e- a0T Q e d0 (8.8) 



9) 



for some c > depending only on a and a. Using Lemma [3.21 we get that as n — > oo 
j e -*e T Q° d o = (i + o (exp(-cn 2£ ))) / e'^^dO. 

Proof of Lemma \6. 6 A Note that 

||P(0)||oc = \\0- 9[1, 1, . . . , 1] T |U < 2||6>|U, (8.10) 

where 

g= gi + ft + ..A (8n) 
n 

Thus 

t4(V 1/2+£ ) c {0 I P(0) G ^(n- 1 / 2 ^)} (8.12) 
Since Q0 = QP(0), using (E2D and <IQ2jl . we get that 



/ exp(-i £ A^W / 



Lnl/„(n- 1 /2+ £) \ / LnC/„(n- 1 /2+ e ) 



| e-^ e d0/ j e-^ 2 dx>^L j 



U3) 



e 2- Q°dO. 

P(0)eC/„(n- 1 /2+ £ ) -oo 

Combining (Q, ([HID and we obtain (l6T25|) . 

To prove Lemma 15.21 we separate the integrand in ( 15 .32 j) into three factors. 
• exp (i T Qcx. n ) — the oscillatory factor, 
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exp J a A 2 k I — the regular factor, 

(vj,v k )eEG 

• exp J b A^ k + R n {0) J — the residual factor. 

\ (v jt v k )eEG J 

The proof consists of the following steps. 

1. In Lemma [8.31 we estimate an integral analogous to ( 15 .32 j) but without an oscillatory 
factor. 

2. Using Lemma [8. 5[ we get rid of the oscillatory factor in ( I5.32p . 

3. Combining Lemma [8.51 and Lemma [8.3[ we complete the proof of Lemma [5.21 
At first, we prove two technical statements. 

Lemma 8.1. For any a > and sequence of functions r n (x) such that as n — > oo 

sup \r n (x)\=o(l). (8.14) 

\ x \< n -l/2+e 

Then as n — >■ oo 

n -l/2+e n -l/2+« 

J nx 2 e- anx2+r " (x) dx = (^ + o(l)^j J e - anx2+rJx) dx. (8.15) 

_ n -l/2+e _ n _l/2+ e 

and 



l/2+e re" 1 / 2 
2 



n 2 x 4 e -anx*+r n (x) dx = J ±_ + Q ^ j j e -anx*+r n (x) fa (g lg ) 



_ n -l/2+e _ n -l/2+e 

Proof of Lemma \8.1\ Note that 

n E +oo n e 

j t 2 e~ at2 dt =(1 + (exp(-cn 2£ ))) J t 2 e~ at2 dt = (J- + o(l)^ j 



e- at dt (8.17) 



and 

y t 4 e~ at2 dt =(1 + (exp(-cn 2£ ))) y tV^d* = (J^ + o(l)\ J e~ at% dt. (8.18) 



Using (j8.14p . we get that as n — > oo 



»■»(*)_ 1| =o(l). (8.19) 



sup \e 

|a.|< n -i/2+ 

Making the substitution t = y/nx and combining ( I8.17P and ( I8.18P with f)8.19p . we obtain 
(I8.15P and (I8.16p . respectively. ■ 
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Lemma 8.2. Under the assumptions of Lemma \8.1\ letr n (x) be differentiable and as n — > oo 

sup |7^(a;)| = 0(n- 1 / 2+3e ). (8.20) 

\x\<n-V 2 + c 

Then as n — >■ oo 

y xe- anx2+Tn{x) dx = O (n- 3/2+4e ) J e~ anx2+rn{x) dx. (8.21) 

_ n -l/2+e _„-l/2+ e 

Proof of Lemma \8.2[ Note that 

n e +oo 

/ te^d* = (1 + O (exp(-cn 2e ))) / te'^dt = 7^ + (exp(-cn 2e )) . (8.22) 
o o 

According to the Mean Value Theorem, we have that for some \x\ < \x\ 

| e rn(«) _ e r n (-«B)| = \ e r ^) r ' n (x)2x\. (8.23) 
Using (I8.20p . we get that as n — > oo 

sup |e r " (x) - e^ - ^ | = O (n- 1+4e ) . (8.24) 

bKn-i/2+s 



We have that 



-1/2+e „-l/2+e 



-anx 2 +r. ' - ' - 



{x) dx = / xe"^ 2 (e rn(x) - e rn( - x) ) dx. (8.25) 



_ n -l/2+e 

Making the substitution t = y/nx and combining (18. 22ft with ( j8.24j) . we obtain (I8.2ip . ■ 
We use notation 

n 

According to the Generalized Mean Inequality, we have that 

Hx/n< Wn) 1/4 . (8.27) 

Since 

fc = (4 + !)#* + e v (8.28) 

and (see lEEIBjl ) 

4 > cm/ 2 (8.29) 

we obtain that 

\0 k \ < — { \4>k\ + } '\0j\ I (8.30) 



an 
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Using Lemma [3. 2\ we find that 



j2m<\m\i<-M\i = -vi 

*— d n n 

3T* 



5.31) 



Combining f[8T2?j) . (1O0]) and flOTj) . we get that 

16 

- ! I-/, I + 

n 



a' n 1 



k\ + —fil) < -t-j (\<f>k\ +Coo(/i 4 /^) 1/4 



a 4 n 4 



Using the inequality (x + y) 4 < 8(x 4 + y 4 ), we obtain that 



fl 4 < r ,Th. + r ^1 
n 4 n b 



.32) 



133) 



where constants c^, c M > depend only on er. 



Lemma 8.3. Let assumptions of Lemma \3.S\ hold. Let {a n } be sequence of positive numbers 
having limit a > 0. Then for any b > as n — >■ oo 



exp ( -a n T Q0 - 6 ]T A 4 fc ] d0 = Q klM ( e - a " eTQf? r/fi 

(vj,v k )eEG 



where constants ki, k% > depend only on a, b and a. 

Proof of Lemma \8. 31 Using the inequality (x + y) 4 < 8(x 4 + y 4 ), we find that 

n 

{v h v k )tEG j=l 

n 

We define R\{B) = 8n^2 9 4 . Thus we have that 



.34) 



.35) 



exp ( -a n T Q0 -b ^ A%]dO> j e^ 6 Q0 - R ^dO. 

(vj,v k )eEG 



Using ( 18. 91) . we find that as n oo 

C/„(n-V2+e) 



+0 (exp(-cn 2£ )) / e- aneT ^ 9 de 
for some c > depending only on a and a. It follows that 



-1/2H 



-1/2H 



/ 



4„-l/2+ e 



-1/2+s _ n -l/2 + e 



C/n(n-i/2+ £ ) 



\ 



\-± n -l/2+e 

</]R 



• • • dd n 



e -*n0TQ9 d0 



.36) 



.37) 



.38) 
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Using (18. 6p . we find that 



-l/2+e 



-1/2+e 



/ ± n -l/2+6 



_ n -l/2+e -„-l/24 
n -l/2+e n -l/2+e 



J ... J J ^-^qe-Rm^ 



\ 



\ 1 i/2+ e 

/ | n -l/2+. 



(i^2 • • ■ d9 r , 



-angi(e 2 ,...,e n )-R 2 (e) 



-1/24 



-1/2+e 



,2 



d0i 



\_4 n -l/2 +e 
\ cr 



(I82 • • • G^r; 



(8.39) 



where i?2(#) = 8n 6^. Using (I8.30p . we get that as n — > oo 



i=2 



4„-l/2+ e 



y ^e-^A-^dO^ (l + 0(exp(-m 2£ ))) ^ (8.40) 

_l n -i/2+ £ |0i|<n 1 /2+ £ 

Combining (I8.38p . ( I8.39p . f )8.40p and Lemma I8~T1 with x = (fii/n, we obtain that as n — >• oo 

g -a„0 T Q0- J R 1 (0) rf + 



/ ^e^" ^-^)^ < c'n 2 [ 

i[/„(n- 1 /2+ e ) if; 



a„(n-V2+e) 

+0 (exp(-cn 2£ )) I e- an0l Qe dO 



for some constants c, c' > depending only on a and a 

Combining similar to f)8.4ip inequalities for <pi,4>2 
as n — )• oo 



, . . . , 



Jl/„(n- 1 /2+ £) 



4 e -a n e T Qe-R 1 (6) d Q < 



(c + c M )c' 



rr 



(8.41) 

and using f )8.33p . we find that 

e -a n T Q0-R 1 (0)^ + 



C/„(n-V2+e) 



+0 (exp(-cn 2£ )) / e- an0T ^ e d6 



(8.42) 



for some c > depending only on a and cr. Note that as n — > oo 



-1/2+e l/2+ £ 



/ n -l/2+ E 



-a„ 9 i(e 2 ,...,e n )-R2(f) 



-1/2+e 



-1/2+e 



\_ n -l/2+e 



/ n -V2+6 

/ e - a «5T+r (l-86n^ + 0(n- 2+8e ))^i 



d02 ■ ■ ■ d9 ri 
\ 



(8.43) 



\_ n -l/2-t 

Combining (18.421) and (I8.43[) . we get that as n — > oo 



<i^2 • • • d9f) 



) 



s 

JU„(n 



e -a n e+Q0-R l{ O) de> 1 



( n -l/2+£) 



e -a n 9 T Q0-R 2 (0) d Q + 

C/„(n- 1 /2+s) 

+0 (exp(-cn 2£ )) / e- an ° T Q e dO, 



(8.44) 
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where c depends only on a, b and cr. 
We continue similarly to (I8.44p 



/ 



e -a n e T Qe-R k (0) dg > 1 1 + £ \ / e -a n e T Qe-R k+1 (0) de+ 

+0 (exp(-cn 2£ )) / e - a " 0T ° rf6/, 



where 

n 

i4(0) = 8nj^ 4 . 
Combining all inequalities of (I8.45P for Ri, R2, . . . , Rn, we get that 

Ju n (n-V 2 + e ) V 71 / il/„(n- 1 /2+ £ ) 

+0 (exp(-cn 2£ )) / e'^^dO 
for some c > depending only on a and a. Note also that 

Combining ([HID, fl8^7j) and flH^Hj) . we obtain fl8TMj) . 



145) 



(8.46) 



147) 



(8.48) 



Lemma 8.4. Let the assumptions of Lemma \3.2\ hold. For positive constants a,b,di,d2 let 
sequence of differentiable functions {R n (6)} be such that 

Re{R n {0))<di 6l(} ° 



n 



an 



dforOe Un^n- 1 ' 2 ^ 



dRJO) 



d9 k 



Then as n — >• 00 



/ 



<d 2 n- 1 / 2+3e . 
( 



C/„(n-V2+ e ) 



\ 



e i±<l> k -aO T QO+R n (0) d Q 



+ 



V„(n- 3 / 2 + e ) 

+ 0{n- 1 ' 2+Ae ) I e^^dO 



and 



( 



f e i±<j> k -ae T Q0+Rn(6) d Q _ q 



ki,k 2 



C/„(n-V2+e) 



7* 



e i±<p k -a0 T Q0+R n (0) d Q 



\ (7„(n-3/ 2 +=) 



+ 



/ 



.49) 



.50) 



.51) 



(8.52) 



+0 (exp(-£n 2£ )) / e - ae Qe d0, 
where constants k\, k%, c > depend only on a, b, d\, d 2 and a . 
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Proof of Lemma\8.4\ For our purpose it is convenient to assume that k = 1. Note that 



\e i fr- a * r & +R »W\dB < J e- aeT ® e+d ^dO. (8.53) 

C/„(n-!/2+c) U n (n-V 2 + S ) 



Using ( 18.2ft . we get that as n — > oo 



J | e i>-a» r <5»+fl»W| d0 = O ( 1 ) y e - a0T ^6>. 

C/„(n-V2+e) l/n ( n -l/2+ e) 



Similar to (I8.37p . for m = 1, 2 we find that as n — > oo 



m\i^- r j} 1 -ae T Qe+R n {e) 



C/u(^n- 1 /2+ e )_[/ n ( n -l/2+ S ) 



|d0 



O (exp(-cn 2£ )) / e- aeTCie dd 



for some c > depending only on a, d\ and a. It follows that 



_„-l/2+c 



-1/2+e 



C?^2 • • • G^n 



/ 



C/„(n-V2+s) 



We define 



and 



\_± n -l/2+e 



r 1 {B) = R n {0)-KQ)- 
According to Mean Value Theorem, for 9 6 U^-nT 1 ^ 6 ) we have that 



0i 



0(n 



-l+4e\ 



Using (|8.6|) . we find that 



-1/2+e -1/2+e 



/ 



ir,- 1 /2+ e 



-1/2+e 



-1/2+e 



\ 



-1/2+e 



0"V£4>1 -« 6 T Q6+Rn (0) ^ i 

\_±„-l/2+e 

4 n -l/2+s 



G?$2 • • • G^ri 



-agi(0 2 ,-,e„)+R' n (0) 



4>i 



\ 



^ e i " rpl - a ^ +ri{0) d9 1 



-1/2+e 



_± n -l/2+e 



d6o . . . d6 r , 



(8.54) 



(8.55) 



(8.56) 



(8.57) 
(8.58) 

(8.59) 



(8.60) 
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Using (I8.30p . we get that as n — > oo 

i n -l/2+ £ 



. i n -l/2+ £ 



(l + 0(exp(-cn 2e ))) 



(8.61) 



|0i|<n 1 /24 



Combining (IH36|) . (KTSOj) . (j8~BT]l with m = 2 and (E35D with x = 0i/n, we obtain ([832]) . 
Note that 



= y 0! (l + z^ a + (^ 1+2£ )) e^^ 1 ^ 
Since dr\jdd\ = dR n /d9i, using (18. 2 1 f) with x = (fii/n, we get that 



(8.62) 



(8.63) 



l0l|<" 1/2+E 



Combining flH36|) . (KTEOj) . flH3Tj) with m = 1 and f l832|) . (gUD with 6 = and (KT55j) . we 
obtain ( KT5T|) . ■ 



Lemma 8.5. Let the assumptions of Lemma \3.S\ hold. For positive constants a,b,di,d2 let 
sequence of vectors {f3 n } and sequence of real differentiable functions {R n (6)} be such that 



||A»||oo < b, 

e T Qe 



R n (0)<d v 



n 



an 



dforOe Un^n- 1 ' 2 ^ 



8R n {0) 



36 



<^n- 1 / 2+3£ . 



Then as n — >■ oo 



(7„(n-i/2+ £ ) 



( r X 

V„(n-l/2+ £ ) 

+0(n- 1 / 2+4e ) /" e'^^dO 



(8.64) 
(8.65) 

(8.66) 



(8.67) 



where R n = R n — | /J 2 ^ 2 and constants k±, k 2 depend only on a, b, d\, d 2 and a 
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Proof of Lemma \8. 51 Using f )8.54p . we get that as n — > oo 

e iPle-aO T Qe+R n {0) d Q = 



l/ n (n-l/2+ £ ) 



_ y e i{Ple-Mi)-aO T Q6+R n (e)-\plel d Q + q ^-3/2+3^ /" e ~aO T Q0 d Q + 

(7„(™- 1/2+£ ) E/„(n-Va+e) 

+ y iABie i^ T n e-P^)-a0^Q0 +Rn {9) de _ 
Taking into account (I8.28P and using Lemma 18. 4[ we find that as n — > 00 



(8.68) 



/ 



{PlO-p 1 6 1 )-aO T Q0+R n {9) jq 



t/„(n-V3+e) 



\ 



\C/„(n-3/2+ £ ) 



+ 



/ 



+0(n- 1 / 2+4£ ) y e- a0T< 3 ^ 



e 



fcl,A:2 



e i( ( 9j0-/3iei)-a0 T Q0+R, l (0)-i/3 1 2 e2 rf g 



\t/„(n-3/2+ E ) 



+ 



+0(n- 1 /2+4, ) y e -«« T 0« d0> 
where constants fci, A:2 depend only on a, 6, di, c?2 and cr. 



According to Lemma 13. 2\ we have that 



\\Q-% = \\Q 



-in ^ 00 

00 

n 



Thus as n — >• 00 



q i (0l0-p 1 e 1 )-a0 T Q0+R n (e) d 



[/n(n-V2+ £ ) 



< 



n 



where c > depends only on a, b, di, 0I2 and a. 



+ 0(n- 3 / 2+4£ ) y e- a0T ^d0, 



(8.69) 



(8.70) 



(8.71) 
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Combining (I8.68P and (18 .71 p . we get that as n — > oo 



C/„(n-V2+s) 



n J 



e i(f3ZO-p 1 ei)-aO T Q0+R 7l (O)-±f3pl d Q + 



C/»(n-3/2+«) 



+ 0(n~ 3 / 2+4£ ) Je^&dO, 



(8.72) 



where Ic^l < c/?i < cb. 
We define 



We continue similarly to f !8.72f) 



k 

E 

3=1 



/3j0j )-a0 T Q0+Rg\e) 



U n (n-V 2 +e) 



1 + 



-(fc+1) 



n 



dO 



i\ 019- ^ )-aO T QO+R { n k+1) (0) 



3=1 



d0+ 



C/„(n-3/2+ £) 



+ 0(n- 3 / 2+4£ ) J e- aeT Q e d0, 



(8.73) 



(8.74) 



where |c (fc+1) | < c/3 k < cb. 

Combining all inequalities of ( I8.74p for k = 0,1, . . . ,n — 1, we get that as n — > oo 



e iPlO-aO T QO+R n (0) d Q 



t/„(n-V2+e) 



1 + 



.(1) 



n 



1 + 



,(ra) 



n 



e - a e' Qe+R n (8) dg+ 

(7„(n-3/2+^) 

+ 0(n~ 1 / 2+4£ ) / e- a0T WdO. 



(8.75) 



Since |c^| < be for = 0, 1, . . . , n — 1, using (I8.75p . we obtain (I8.67p . 



Proof of Lemma \5.^ . Note that for G U n (-n 1 / 2+£ ) as n — > oo 



M S A ^ 



0{n 



-l/2+3e\ 



(8.76) 
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We define (3 n = Qct n . Using Lemma [8. 5 [ we find that as n — > oo 

J exp j i T Qa n - aO T Q6 -b A % + R n{0) J dO 

n -i/2+e) \ (vj,v k )€EG J 



U„{ 




ki,k2 



\ 



exp l -aO T QO - b ^ A % + ) d0 

+0(n 



+ 



-l/2+4e^ / e -aO T Q8 



(8.77) 



where R n = R n — | ^2 and constants fei, A)2 depend only on a, b, di, di and a. Note that 



for some 0^3 > 0, depending only on c and a, 



^.75 



Combining f)8.2p , (I8.77P and Lemma 18. 3[ we find that as n — > 00 

J exp j i 6 T Qa n - aO T QO -b A % + R n{Q) J dO 



(8.79) 



n — 1 n — 1 



®k[,k 2 (^7r^/ydet<^ 

where constants k[, k! 2 depend only on a, b, d\, g?2 and a. Note that 

+00 

j e~ anx2 dx j exp I i T Qa n -a £ A? fc - 6 A i + ) rfL 



Lnvb 



(vj,v k )eEG 



(vj,v k )eEG 



J exp j z T Qa n - a0 T Q0 - 6 Aj fc + i? n (0) j dO. 



(8.80) 



p(6»)eLnv 



We have that 



exp I 1 6 T Qa n - a0 T Q0 - b A % + Rn(P) 

{vj,v k )eEG 



< e -ae T Q°+^o T Q0. (8.8i) 



Thus, combining ([HID, (1532} , (157731) and ( 18780]) . we obtain fl5733|) 
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9. Final remarks 

In fact, using Lemma [5. II and Lemma 16.31 the estimation of the number of Eulerian circuits 
is reduced (see proof of Lemma 15. 3ft to estimating the integral 

e W Ue T Q<x-\ £ £ A% + ±tr(AQ- l )AdO, (9.1) 

\ (vj,v k )eEG (vj,v k )£EG J 

where a denotes the vector composed of the diagonal elements of Q~ l , A denotes the diagonal 
matrix whose diagonal elements are equal to components of the vector Q6. Apparently, it is 
possible to estimate integral (19.1 j) more accurately for particular classes of graphs and obtain 
asymptotic formulas for Eul(G), similar to ()2.12p . 

Finally, we want to note that the following expression 

2^-^n-^^mf[(^-l)\ (9.2) 

j= i V / 

gives a surprisingly good estimate for the number of Eulerian circuits in graphs. Namely, we 
calculated the exact numbers of Eulerian circuits for small random graphs and in all cases 
the values given by (19.21) differ from the exact ones within not more than 30% error. 
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